Renormalisation Flow and Universality for Ultracold Fermionic Atoms 
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A functional renormalisation group study for the BEC-BCS crossover for ultracold gases of 
fermionic atoms is presented. We discuss the fixed point which is at the origin of universality 
for broad Feshbach resonances. All macroscopic quantities depend only on one relevant parameter, 
the concentration akp, besides their dependence on the temperature in units of the Fermi energy. 
In particular, we compute the universal ratio between molecular and atomic scattering length in 
vacuum. We also present an estimate to which level of accuracy universality holds for gases of Li 
and K atoms. 
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I. INTRODUCTION 

The quantitatively precise understanding of the 
crossover from a Bose-Einstcin condensate (BEC) to BCS 
superfluidity [l| in gases of ultracold fermionic atoms 
2, 3] is a theoretical challenge. If it can be met, the com- 
parison with future experimental precision results could 
set a new milestone for the understanding of the tran- 
sition from known microscopic laws to macroscopic ob- 
servations at length scales several orders of magnitude 
larger than the characteristic atomic or molecular length 
scales. Furthermore, these systems could become a ma- 
jor testing ground for theoretical methods dealing with 
the fluctuation problem in complex many-body systems 
in a context where no small couplings are available. The 
theoretical progress could go far beyond the understand- 
ing of critical exponents, amplitudes and the equation of 
state near a second-order phase transition. Then, a whole 
range in temperature and coupling constants could be- 
come accessible to precise calculations and experimental 
tests 0. 

The qualitative features of the BEC-BCS crossover 
through a Feshbach resonance can already be well re- 
produced by extensions of mean field theory which ac- 
count for the contribution to the density from di-atom or 
molecule collective or bound states [||. Furthermore, in 
the limit of a narrow Feshbach resonance the crossover 
problem can be solved exactly [|| , with the possibility of a 
perturbative expansion for a small Feshbach or Yukawa 
coupling. A systematic expansion beyond the case of 
small Yukawa couplings has been performed as a 1/N- 
expansion Q- The crossover regime is also accessible 
to e-expansion techniques @. At zero temperature Q 
as well as at finite temperature (Toj . numerical results 
based on various Monte-Carlo methods are available at 
the crossover. A unified picture of the whole phase di- 
agram has arisen from functional field-theoretical tech- 
niques, in particular from self-consistent or i-matrix ap- 
proaches fill . Dyson-Schwinger equations @, 2PI 
methods [13j, and the functional renormalisation group 
(RG) M- 

Recently, it has been advocated [6| that a large uni- 



versality holds for broad Feshbach resonances: all purely 
"macroscopic" quantities can be expressed in terms of 
only two parameters, namely the concentration c = akp 
and the temperature in units of the Fermi energy, T — 
T/ep. Here, kp and ep are defined by the density of 
atoms, n = fc F /(37r 2 ), ep = fc|/(2M) and a is the scat- 
tering length. Universality means that the thermody- 
namic quantities and the correlation functions can be 
computed independently of the particular realizations of 
microscopic physics, as for example the Feshbach reso- 
nances in 6 Li or 40 K. A special point in this large uni- 
versality region is the location of the resonance, c — * oo, 
where the scaling argument by Ho applies. 

We emphasise that the universal description in terms 
of two parameters holds even in situations where the un- 
derstanding of the microscopic physics may necessitate 
several other parameters. It is thus much more than the 
simple observation that the appropriate model involves 
only two parameters. The fact that two effective param- 
eters are sufficient is related to the existence of a fixed 
point in the renormalisation group flow; for a given T, 
this fixed point has only one relevant direction, namely 
the concentration c. Then, the question of accuracy of 
a two-parameter description is related to the rate how 
fast this fixed point is approached by the flow towards 
the infrared. It depends on the scales involved and the 
physical questions under investigation. As in statistical 
physics, the deviations from universality can be related to 
the scaling dimensions of operators evaluated at infrared 
fixed points. We will address these questions quantita- 
tively for Feshbach resonances in 6 Li or 40 K. 

Even for broad Feshbach resonances not all quantities 
admit a universal description. An example is given by the 
number density of microscopic molecules which depends 
on further parameters 0, 03, Gil • 

The quantitative study of this universality has first 
been performed through the solution of Schwinger-Dyson 
equations, including the molecule fluctuations (a. [l2l \vi\ . 
Already in this work it has been argued that the key for a 
proper understanding of universality lies in the renormal- 
isation flow and its (partial) fixed points. A similar point 
of view has been advocated by Nikolic and Sachdev [7[. 
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A first functional renormalisation group study evaluat- 
ing the renormalisation flow for the whole phase diagram 
has been put forward in [3]. In the present paper, we 
perform a systematic study of the universality aspects 
of the renormalisation flow for various couplings thereby 
extending the results of [3] concerning universality. We 
relate the universality for broad Feshbach resonances to 
a non-perturbative fixed point which is infrared stable 
except for one relevant parameter corresponding to the 
concentration c. 

Our approach is based on approximate solutions of ex- 
act functional renormalisation group equations for 
reviews see 0, [2(| [Hj], an d f° r applications to non- 
relavtivistic fermions sec [26]. These are derived by vary- 
ing an effective "infrared" cutoff associated to a momen- 
tum scale k (in units of fcp). The solution to the fluc- 
tuation problem corresponds then to the limit k — > 0. 
In this way, a continuous interpolation from the micro- 
scopic Hamiltonian or classical action to the macroscopic 
observables described by the effective action is achieved. 
In this paper, we will use a very simple cutoff associated 
to an additional negative chemical potential. This works 
well on the BEC side of the crossover and in the vac- 
uum limit of vanishing density and temperature. Since 
the crucial characteristics of the fixed points relevant for 
universality can be found in the vacuum limit, this will be 
sufficient for our purp ose. The investigation of the whole 
phase diagram in [I4| has been performed with a different 
cutoff more amiable to this global task, see [2l|, [22j • 

The simple form of the cutoff allows for analytic solu- 
tions of the flow equations; these can be used for a di- 
rect check of the more general arguments, resulting from 
the investigation of fixed points and their stability. As 
a concrete example, we compute the universal ratio of 
the molecular scattering length to the atomic scattering 
length in vacuum. This ratio is reduced as compared to 
its mean field value by the presence of fluctuations of 
collective bosonic degrees of freedom. We also present 
quantitative estimates to what precision universality is 
realized for ultracold gases of Li and K. 

In Sect. [Til we introduce the functional-integral formu- 
lation of our model for the Feshbach resonance. It con- 
tains explicit bosonic fields for di-atom states [23l.l24l.l25l]. 
In the limit of a divergent Feshbach or Yukawa coupling 
h v , this model is equivalent to a purely fermionic formu- 
lation with a point-like interaction. In Sect. lIIll we briefly 
recapitulate the framework of the exact flow equation for 
the effective average action, and specify our cutoff and 
truncations. The initial values of the flow at the micro- 
scopic scale are given in Sect. HVl while explicit formulae 
for the flow of couplings and the effective potential are 
computed in Sect. IVl and I VII 

In Sect. IVII1 we specialise to T = and solve the flow 
equations for various couplings explicitly within our trun- 
cation. Sect. IVIIII generalises these results and presents 
a general discussion of the fixed points in the flow of the 
rescaled dimensionless parameters /i^,, and A,/,, corre- 
sponding to the Yukawa coupling, the energy difference 



between open- and closed-channel states and an addi- 
tional point-like four-fermion vertex which accounts for 
a "background scattering length" . In particular, the fixed 
point describing broad Feshbach resonances has only one 
relevant parameter c, whereas the fixed point characteris- 
tic for narrow Feshbach resonances has two relevant pa- 
rameters c and hfp. Sect. IIXI supplements a discussion 
of the density which is needed for a practical contact to 
experiment, and Sect. IXl relates the initial values of the 
flow (the classical action) to observable quantities, such 
as the molecular binding energy or the scattering length 
depending on a magnetic field B. 

In Sect. IXIl we discuss the fixed-point behaviour for 
an additional parameter, namely the effective four-boson 
vertex X v . This fixed point is responsible for the uni- 
versal ratio between the scattering length for molecules 
and for atoms. We finally present a general discussion 
of universality for ultracold fermionic atoms in Sect. lXlTl 
and we estimate the deviations from the universal broad- 
resonance limit for Li and K in Sect. IXIIII Sect. IXIVI 
contains our conclusions. 



II. FUNCTIONAL INTEGRAL 

We start from the partition function in presence of 
sources (J dx — J d 3 xdr) 

ZbM= fwexp{-S[^^] (I) 

+ / <®[jfa)to) + Mx)4* (*)]}, 

with 



S = dx 



A 

+<f>* (d T - — + v + bv - 2a) <t> 



(h v + Sh v ) ( (t>*ipiip2 — <ptp*^2 



Here, we have rescaled the fields and couplings, to- 
gether with the space and time coordinates x = kx, 
f = (k 2 /2M)t, with M being the mass of the atoms. Our 
units are h = c = ks = I. We use the Matsubara formal- 
ism with Euclidean time r on a torus with a circumfer- 
ence given by the inverse temperature T . The thermo- 
dynamic variables are T — 2MT/k 2 and a — 2Ma/k 2 , 
with a denoting the effective chemical potential. 
The model parameters are: 

(1) the detuning of the magnetic field B — Bo (with 
fl = 2fi B for 6 Li and p = 1.57/i S for 40 K) 



2M , 
v=-^-p(B-B ) 



(2) 
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(2) the Feshbach or Yukawa coupling h v which can be 
extracted from the properties of the quantum me- 
chanical two-atom system as the molecular binding 
energy or the scattering cross section, 

(3) a point-like interaction for the fermionic atoms pa- 
rameterised by X^. The shifts Si>, SX^, and 5h<p are 
counter terms that are removed by the ultraviolet 
renormalisation. Details of the rescaling and the 
formulation can be found in [l^, [H| • 

The scale k (which sets the units) is arbitrary. (A typical 
value is k = leV.) Under a rescaling of k — > all 
quantities scale according to their canonical dimension, 
i.e., 



x 
f 



X 



= kx — > xj \x 
-> [i 2 f, a- 



(k 2 /(2M))T^r/^, 



/x 2 i>, 



3/2, 



(3) 



We observe that the canonical dimension of time is mi- 
nus two and therefore the nonrelativistic Lagrangian has 
dimension five and not four, as for a relativistic quantum 
field theory. The total atom density n defines the Fermi 
momentum fcp, 



fcp 



3tt 



2 ' 



(4) 



If one associates k with kp all quan tities are expressed in 
units of the Fermi momentum [l2| . 

From the partition function the effective action F ob- 
tains by the usual Legendre transform, (p = ((f>), 



(5) 



T[ip] = -InZU^ip]] + / dx(j%(p + j^tp*) 



The order parameter tpg for superfluidity corresponds to 
the minimum of Y for = and obeys the field equation 



ST 

60 



= 0. 



■pa 



(6) 



The effective action generates the 1PI Green's functions 
such that the propagators and transition amplitudes can 
be directly related to the functional derivatives of T. The 
construction above is easily extended to an effective ac- 
tion that is also a functional of fermionic fields by intro- 
ducing the appropriate fermionic sources, see below. 



III. EXACT FUNCTIONAL FLOW EQUATION 
AND TRUNCATION 

The variation of T with the change of an effective in- 
frared cutoff k is given by an exact renormalisation group 
equation [IH, Q1J [2(| HH • For the present theory, this ap- 
proach has been implemented in [l4| within a study of 
the phase diagram; see also 26]. For RG studies of the 
purely bosonic system, see [271 ] . 



For the purpose of the present work, we extend the 
truncation used in [3] as well as using a particular ver- 
sion of the flow equation where the cutoff acts like a shift 
in the respective chemical potentials for ip and (p, as also 
done in [281 ]. To that end, we introduce an infrared- 
regularised partition function Z k by adding a cutoff term 
to the action in Eq. fTJ), 



S$,$^$[tM]+AS fc hM], 



where the cutoff term AS k is chosen as 



AS k [ij,4>] = I dx{R { pi>H + R k a> (t> 



(7) 



(8) 



with 



R { P = Z^(k)k l , Rf = 2Z v (k)k\ 



OR 



(F) 



dk 2 



QlpZ^J; , 



OR 



(B) 



dk 2 



d 



(9) 



The fc-dependent wave function renormalisation Z^, ttf) will 
be determined below, and we note that k is measured in 
units of the fixed scale k. Introducing the effective aver- 
age action in terms of a modified Legendre transform, 

r fe $, 0\ = - In Z k [ty , j^] - ASk [i>, 0) 

+ J dx(j;0 + 3<P ip* + - ft fa), (io) 

the exact functional renormalisation group equation (flow 
equation) can straightforwardly be derived, 



dk 2 



dx{Z^Q^{ft%li) c + 2Z V Q V> (( I 
Tr{-Z^Q4T k 2) +R k )^ 
+2Z v Q v (T ( k 2) +R k )- 1 ,.}. 



(11) 



Here, we have expressed the propagators of the fermionic 
and the bosonic fields by the corresponding components 
of the inverse of the matrix of second functional deriva- 
tives of T k . The trace Tr contains an integration over 
x or a corresponding momentum integration as well as a 

(2) 

trace over all internal indices. Both T^. and are func- 
tionals of arbitrary fields (p and ip which are kept fixed 
for the k derivative in Eq. (fTTj) . 

The flow equation (jlip is a nonlinear functional differ- 
ential equation and we can only hope to find approximate 
solutions by suitable truncations of the most general form 
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of Tk- In this work, we exploit the ansatz 



d.r 



+Z v ip* (d T — A v A) + u(0, a) 



= dx 



+f*(d T - A v A)ip + u(ip, a) 



(12) 



In the second equation, we have used renormalised fields 



,1/2 



h v = K /2 



and renormalised couplings 



hip — h v Z^ 1 Z ip 1 / 2 , 



— A^Z^ 2 . 



(13) 



(14) 



The ansatz (TT2"j) generalises that of [l4[ by the four- 
fermion interaction with coupling A.0. This term is par- 
ticularly important in the limit of broad Feshbach reso- 
nances as will be discussed later. 

The global symmetry of U(l) phase rotations implies 
that the effective potential u can only depend on p = 
f*tp. If the ground state corresponds to a homogeneous 
tpo the pressure and density can be computed from the 
properties of the effective potential 



du 

dip 



, uo(o-) = u((po,&) , 



P = 



' n = -P dU ° 



2M 



da 



(15) 



At this point, some comments concerning the wave 
function renormalisations are in order: 

(i) Analytic continuation to "real time" and a Fourier 
transform to (real) frequency space results in d T — > 
—lo. We define Z$ tV by the coefficient of the term 
linear in to in the full inverse propagator P\j,,tp given 
by the terms quadratic in ip or in T^. More pre- 
cisely, we choose in Fourier space 



Z v — 



q=0,U)=0 



(16) 



(ii) With this definition of Z^ and Z V1 the renormalised 
fields tp, if have a unit residuum for the pole in the 
propagator for q — * if the "on shell value" of oj 
vanishes for q = 0. 

(iii) We use the same Z^ m in the definition of the cutoff 
© as in Eq. Q2J). 



(iv) The fields ip and describe "microscopic" or "bare" 
atoms and molecules while the renormalised field 
(f describes dressed molecules. The wave func- 
tion renormalisation Z v accounts for a description 
of dressed molecules as a mixing of microscopic 
molecules and di-atom states 0, 0, [13] • 



IV. INITIAL CONDITIONS 



Inserting the truncation (|12p into the flow equation 
(fTT) . and taking appropriate functional derivatives, leads 
to a coupled set of differential equations for the couplings 
Z^jA^jZ^jAipjhipjXt/, as well as the effective potential 
u(ip). The task is to follow the flow of these couplings as 
the cutoff scale k is changed. The initial values for large 
k are taken as 



Z(p — Zip 



1 



Aip — — , Aty — 1, 



m 



u = m v p 

2 



v + bv — 2(T. 



(17) 



Here v is related to the detuning of the magnetic field 

(18) 



and we have to fix 8v (i.e., the renormalised counterpart 
of 8v in Eq. (TTJ) such that the Feshbach resonance in 
vacuum occurs for v = 0. The initial values of and h v 
are chosen such that the molecular binding energy and 
scattering of two atoms in vacuum is correctly described. 

A realistic model contains an effective ultraviolet cutoff 
A which is given roughly by the inverse of the range of 
the van der Waals force. For k -C A/fc we can take the 
limit A — > oo since all momentum integrals are ultraviolet 
finite. All ultraviolet "divergencies" are already absorbed 
in the computation of Ta |12l . fl4| . For 6 Li or 40 K and 
k = leV one has A/fc of order 10 2 to 10 3 . 



V. RUNNING COUPLINGS AND ANOMALOUS 
DIMENSIONS 

The flow equations for the various couplings and the 
densities derived from (jlip have a simple interpretation 
as renormalisation group improved one-loop equations 
[TH with full propagators and dressed vertices, but are 
exact. The renormalisation constants Z v , Z^ are related 
to the dependence of the unrenormalised inverse propa- 
gators (for and ip) on the (Minkowski) frequency, i.e., 
the coefficients of the terms linear in uj(= — d/dr) in 
r^ 2 ). The second derivative of Eq. (fTTI) y ields exact flow 
equations for the inverse propagator [18L liH l20l . [2lj and 
therefore for Z^jZ^p. We define 



fjip,(p = rj^ iV /2k 2 = -d\nZ^,^/dk 2 



(19) 
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and obtain in our approximation 



hlQ, h d 



2 dfc 2 v ' v 
Here the fermion loop integral 



(20) 



1 f d 3 q 

= ^ J ^3 77^ 3 [tanh 7v - 7(/3 cosh~ 2 7(p ] (21) 
involves 



1 



7 = 4(^<f + fc 2 -a). 



(22) 



In our truncation, the boson fluctuations contribute to 
fj v only in the superfluid phase (fj^, ' (p = 0) = 0) [141 ]. 

Crucial quantities for the investigation of this problem 
are the running of the Yukawa coupling h v , 

= {rj v + 2f^)h% + 2h%\^Q^ F \ (23) 



with 



■/" 1 = I ^77; 5 [(3 7 2 ~7^tanh 7 , 



16T 2 7 (2tt 



+ 7v cosh 2 7 ^(- 3 7 2 + 7 2 

+2 7v ( 7 2 - 7 2 )tanh 7v )], (24) 

and the point-like four-fermion interaction 

A A v , = 2fa + A 2 (If) + J A J , (25) 



4T 2 ./ (2?r) 



3 77y 1 tanh 7v cosh 2 7v . 



For both quantities, we neglect the contribution from 
the molecule fluctuations. This is a valid approximation 
for the parameter ranges for which we give quantitative 
results. In other ranges, the molecule fluctuations give 
more dominant contributions, see [l4| and the discussion 

(-F) (F) 

below. For p — 0, we note Jtp —lip . 



VI. FLOW OF THE EFFECTIVE POTENTIAL 

Keeping now p (and not p) fixed, the flow of the effec- 
tive potential obeys, 

d 

^u(p) = (F(f) + (B(p) + t} v fn*'(p), (26) 



with 



^ I/O = -(),■ / ^^(^tanh^-l). 



Cb(p) = Q, 



d 3 q { a + K2 



(2tt) 3 V a 



cotha^ — 1 



). (27) 



Here, we define 



1 

2T 



{A v q 2 + 2k 2 +> 



{A v q + 2k l 



2pu") 



1/2 



a 



K-2 



A^q 2 



2k 2 



pu 
2T 



2T 

K3 = 



pV" 
2T 



(28) 



and the primes denote derivatives with respect to p. The 
contribution of the boson fluctuations is computed in a 
basis where the complex field ip is written in terms of two 
real fields with inverse bosonic propagator 



A v f 



2k 2 + u' 
2itnf 



2pu" 



-2-kuT 



A v q z + 2k' 

Then, Eq. (fTT|) evaluated for constant <p becomes 
d 3 q _ 1 



C B (p)^fQ v tv^2 J 



(2tt)3 



A v q 2 



2k 2 



(29) 



(30) 



in accordance with [29(. The subtraction of a p- 
independent part renders the momentum integral ultra- 
violet finite. For extensive numerical studies one would 
rather rely on cutoff choices as in [14( which do not re- 
quire any ultraviolet subtraction, and minimise the nu- 
merical costs. However, this is of no importance for the 
quantities computed in this paper - the dependence of A v 
on a and T is negligible, and simplicity of the approach 
is the more important property. 

In the symmetric phase (SYM), one has po — 0, u'(0) = 
m 2 ,u"(0) = \ v , whereas in the presence of sponta- 
neous symmetry breaking (SSB) we use po > 0, u'(po) = 
0,u"(po) = X v . Eq. ([26)1 yields in the symmetric phase 



A TO 2 = h 2 v Q^(p = 0) - Q^fHp = 0) + f) v m 



Here, we define the boson loop integral 



t(b) 




d 3 c 



(2tt)3 

(2a + n 2 )u" + apu 1 ' 
k 2 (k2 - a) - an 3 



(31) 



(32) 



a - 



K 2 



a 



, sinh 



coth a v 



a 



such that 



4^ = 0) 



2jp 
f 



d 3 q . 2 
sinn a. 



(27T) 



(33) 



The flow starts with a positive to 2 for large k and is 
therefore in the symmetric regime. As k decreases, m 2 
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decreases according to Eq. (|31l) . with fj v being positive. 
At high temperature, m 2 , stays positive for all k and the 
system is in the symmetric phase and ungapped for all 
momenta. At temperatures below a pseudo-critical tem- 
perature T p , m~. hits zero at some critical k c . For k < k c , 
the flow first continues in the broken-symmetry regime 
with nonzero pq. Whereas fermionic fluctuations tend 
to increase po, molecule fluctuations cause its depletion. 
For temperatures inbetween the critical temperature and 
the pseudo-critical temperature, T c < T < T p , molecule 
fluctuations eventually win out over the fermionic fluc- 
tuations and po vanishes again at smaller fc; here, the 
system is in the symmetric phase. Identifying the k de- 
pendence of the fermionic 2-point function with the mo- 
mentum dependence, the non-zero value of po for finite 
k can be associated with a pseudo-gap [30(. Below the 
critical temperature T Cl po stays finite for k — > 0, cor- 
responding to the superfluid phase with a truly gapped 
spectrum. 

Molecule fluctuations are particularly important dur- 
ing those stages of the flow where po is non-zero, i.e., 
for k < k c and below the pseudo-critical temperature T p , 



see [14| . In this work, we focus on the properties of the 
flow in the symmetric regime, where po = for all val- 
ues of k under consideration. We still include molecule 
fluctuations for all quantities that are associated with 
the effective potential u(p), in particular the molecule 
density and the molecule-molecule scattering length, see 
Sects. IIXI and IXII below. But molecule fluctuations are 
neglected for the other running couplings which are dom- 
inated by fermion fluctuations in the parameter ranges 
considered here. 



VII. SOLVING THE FLOW EQUATION FOR 
ZERO TEMPERATURE 

Let us first concentrate on T = and k 2 — a > 0. We 
define 

K 2 = k 2 -a (34) 

and employ d/dk 2 = d/dK 2 . Then the flow equation for 
u takes the explicit form 



dK 



:u{p) = fj^pu 



I 

2^ 



oo 




/ dxx 2 I 




1 





A^x 2 + K 2 



-Qu 



yJ(Ai,x*+K*)* + h*p 
( A v x 2 + 2k 2 



(35) 



pu 



^(A^x 2 + 2k 2 + u')(A v x 2 + 2k 2 +u' + 2pu") 



- I 



For p = 0, T = 0, <t < 0, one also finds 



The solution of the last equation, 



t(f) _ J_a-W k -i 



h„ = 0, 



(36) 



and, in our approximation, the contributions of the boson 
loops to fj v , dh^/dk 2 and dX^/dk 2 vanish for m 2 , > 0. 

At T = and a < 0, we can use — A^ — I, = 1. 
This is an exact result [3] as long as all propagators in 
the relevant diagrams have simple poles in the imaginary 
qo plane. With 



64tt 



(37) 



we find the coupled system of flow equations for the sym- 
metric phase 



dK 

dX-4, 
dK 



32?r 
8tF' 



4tt 



(38) 



~~ A i>,o 



+ — (K m - K) 



1 

K 
8^' 



(39) 



renders A^ almost independent of K if K <§; K ln + 
8ir/\^ z - m . Here, we assume implicitly that X^ t - m is not too 
much negative such that X^ remains finite in the whole 
k range of interest. For positive A^,, the self-consistency 
of the flow requires an upper bound A^ !] o/(87r) < l/K- ln . 
In contrast, for A.0 = the solution for h 2 , 



is dominated by small K. 

For X^p ^ 0, the flow of is modified, however, with- 
out changing the characteristic behaviour for if — > 0. 
Indeed, in the limit K — > the term ~ X^ becomes sub- 
dominant for the evolution of h^. The flow for the ratio 
h^/K reaches a fixed point 32-7T. We can explicitcly solve 
the system (|3"5|) by consecutive integrations. One obtains 
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for h 2 v 



Z v h% 



K = 



exp ■ 



1 

47T 



if in 



dxXrf (x) 



K 



h 2 



i 2 9 . m (1 - c K in ) z (1 - c K) 

This yields for the wave function renormalisation 

k- 1 

Z v = 1 + ^ f dix-^hlix- 1 ) 



with 



h 2 
32^ 



-2cq In 



IfK - c 
1 

K 



co 



(42) 



The counterterm <5i> is determined from the condition 
m^B = B ,a = 0, k = 0) = 0, such that (up to minor 
corrections) 



2Mfl(B - B a ) 



h 2 



K 



8?r 1-cqA 



- 2<7. 



(49) 



(41) In particular, for k = and <r < one has K 



and Eq. (|49p yields the expression for rh 2 in the presence 



of all fluctuations 



2M}i(B - B ) 
k 2 



8tt 



- 2a. 



(50) 



At this point, we have the explicit solution of the flow 
equation for T = 0. In the present simple truncation, we 
expect that the description of the system is qualitatively 
correct. Higher quantitative precision will require a more 
extended truncation; however, this is not the main em- 
phasis of the present work which rather concentrates on 
the structural properties related to the fixed points. 



co 



K 



> Cq. 



(43) 



For K 



0, the ratio h^/h 2 - m obviously reaches a con- 



stant which depends on X^^ n , K ln , whereas Z v diverges 
~ h% /(327rir), as is consistent with the fixed-point be- 
haviour h 2 v m 327T. Neglecting Co (as always appropriate 
for small enough K) and assuming a broad resonance 
h 2 _ ^> 327r, we arrive at the simple relation 



h 2 

3271-iT 



(44) 



In this limit, the factor Q v appearing in the flow equa- 
tions is given by 



o» = 1- 



2(k 2 



(45) 



and approaches 1/2 for k 2 ^> —a. We note that in this 
range the bosonic cutoff function Rfjf = 2Z lp k 2 is effec- 
tively linear in k and not quadratic. 

We finally investigate the flow equation for m^, in the 
symmetric phase. As we will see below, the boson loops 
(the last contribution in Eq. (|35|) ) vanish in our trunca- 
tion. This yields 



di 



h 2 

<p 



dK 8tt 



32tt 



-K- 



or 



dm 2 v _ 
dK 

The solution reads 
„ 2 2Mp(B - B 



8 

dK 



K 

8n' 



Km 



k 2 
(i 



2a + 50 



(46) 



(47) 



(48) 



c K in ) 2 



inco 



1 



1 



1 - Coif in 1 — CqK 



VIII. FLOW EQUATIONS AND FIXED POINTS 

The explicit solution of the preceding section is indeed 
very useful for verifying and explicitly demonstrating the 
general fixed-point properties. The following features of 
the flow equations hold in a much wider context of differ- 
ent microscopic actions and different cutoffs. The overall 
pattern of the flow is governed by the existence of fixed 
points. Some of these fixed points may correspond to 
particularly simple situations, being less relevant for the 
physics under discussion. By contrast, the stability or 
instabilities of small deviations from the various fixed 
points are much more important, as they determine the 
topology of the flow in the space of coupling constants, as 
demonstrated in Fig.[T] If the system is in the vicinity of 
any of the different fixed points the number of effective 
couplings needed for a description of the macrophysics 
(beyond T) corresponds to the number of relevant direc- 
tions. 

A particularity of our system concerns a certain redun- 
dancy in the description: a pointlike interaction can be 
described either by the four-fermion coupling A^, or by 
a limiting behaviour of the scalar exchange. Indeed, for 
^.in — * 00 an d fixed hpj n /m?- in the molecule exchange 
interaction becomes effectively point-like. Therefore, we 
may define an effective point-like coupling 



V>,eff 



K hi 
Ai/. — ry — A,/; j (51) 



which describes the interaction in the zero-momentum 
limit. Combining the flow equations (|38p. (|46|) for 

A,/,, hy, ml, one obtains 



dK ~ 8tt ' 



(52) 
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This is the same flow equation as for A^ (|38[). 

Next, we may consider the renormalised couplings in 
units of kk instead of k. This will reveal the relevant fixed 
points for the flow more clearly. We define, according to 
the scaling dimensions of Eq. (J3j) , 

Ay = AyA: , Ai/, >e ff = Aip^&k, 

r ""P ~ 2 2 /; 2 

Infe/fein, (53) 
and obtain the dimensionless flow equations 



0.2 


-0.2 
-0.4 
-0.6 
-0.8 
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-1.2 



1 ' T 
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(D) 



^ f V 
/» f V 



v T A 
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Aw, + 



3 ^ 



8tt 

^y 3 , ^v-y 



8tt 



v 32tt 4tt ' 
v 8tt 32tt ' 



FIG. 1: Location of fixed points (A)-(D) projected onto 
the plane which is spanned by the couplings \^,/(8ir) and 
hp/(32ir). The arrows characterise the flow of the couplings 
towards the infrared. 



with 



V = 



k 2 



1/2 



(54) 



(55) 



we can identify 



MirMfi{B - Bo) 

k 2 h 2 h- 
K n tp,0 hj ™ 



(60) 



Let us first consider a = 0, or, more generally, k 2 >• 

-a, such that y = 1. In this case, we observe two fixed UM , U (n)ml ( 4 , „ I1U „ . ;mu , X( . 

points for Ay,, either Ay, = or Ay = -8tt. The corre- fh%~B-B Q the relevant coupling 
sponding fixed points for hz, ^ are 



We conclude that broad Feshbach resonances (large 
enough Yukawa couplings) can be characterised by the 
fixed point (A), with h v and A^ being irrelevant and 



(We 
, 



(A) : 

(B) : 

will 

m£ = 



A,/, = 



A 



y, 



-8tt 



h% = 32tt , rh 2 v = 4 
= 96tt , 



-12. 



(56) 



discuss later the fixed points with 
= 0.) The fixed point (A) is infrared stable 



h 2 



For the second fixed point (B), Ay, becomes a relevant 
parameter. For Y? , the fixed point (|56|) with hi = 96tt 
remains infrared attractive, but the fixed point for mz 
occurs for negative ml where our computation in the 
symmetric phase is no longer valid. 

Finally, we turn to the fixed points with h = 0. In 

this case, hz, is always a relevant coupling and increases 



for Ay, and h 2 v - both couplings run towards their fixed- as - 1S lowered. The fixed point 



point values as k is lowered, as is visualised in Fig. [T] In 
contrast, ml is infrared unstable - the detuning B — Bq 
corresponds to a relevant perturbation of the fixed point. 
With hi at the fixed-point value, ml deviates from its 
fixed point with an anomalous dimension 



(C) 



A,,, = , h?=0, m 2 =0 



(61) 



~ 2 a i r ~ 2 "'in 

m v = 4 + dm in — . 



(57) 



Comparing this with the explicit solution (|4"9"|) for Ay 

0, a = 0, h v = htpfl, namely 



fh,„ 
v 8tt 



h 2 v 2Mji{B — B Q ) 



Z v k 2 k 2 



(58) 



is infrared attractive for Ay. For this case of a narrow 
Feshbach resonance, the crossover problem can be solved 
exactly 0] , and perturbation theory around the exact so- 
lution becomes valid for small hz,, Ay. The coupling mt 
is a second relevant parameter around this fixed point. 
As one increases h^L , a crossover to the fixed point (A) 
occurs |6J. The "narrow-resonance fixed point" (C) de- 
scribes the limit of a combined model with free fermions 
and free bosons, sharing the same chemical potential. 
We emphasise that, for small h v , the equivalent purely 
fermionic description typically has a nonlocal interaction. 
For the fourth fixed point of our system, 



with Eq. dig) , 



(D) 



A,, 



, m 



T,2 



0, 



(62) 



7,2 

'V,o 
327rfc' 



all three parameters are relevant. This point corre- 

(59) 

sponds again to strong attractive interactions between 
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the fermionic atoms. The flow away from this fixed point 
for non- vanishing hz, /rnt can be characterised by the flow 
of the ratio (y = 1) 



ml 



)- 






4- 








4tt j 




8tt 





(63) 



which does not have a fixed point for positive h 2 ,rrq p . 
Formally, the flow runs for negative m 2 towards the fixed 

point (B) with h^/rh 2 , = Sir. Actually, it may be possi- 
ble to consider A,/, as a redundant parameter using partial 
bosonisation and rebosonisation during the flow [3l| , see 
also 211 ] . Then, the fixed points with A^ = — 87r may 
again be associated with broad Feshbach resonances. We 
show the four fixed points A, B, C, D and the infrared 
flow of the couplings h 2 ^ and A v in Fig. Q] 

We should emphasise that the inclusion of the omitted 
contributions from boson fluctuations to the flow of A,/, 
could result in corrections (y = 1), 



\ 2 

C t A^ — + h + -^~^n 

an 6Zn 



32tt 



(64) 



such that the flow of A^, remains no longer independent of 
h 2 ^ . This will change the precise location of fixed points 
(A) and (B). We expect that the qualitative characteris- 
tics of the fixed point (A) remain unchanged, whereas the 
fate of the fixed point (B) is less clear. We may consider 
the flow of the ratio X^/h 2 ,, 



(65) 





c ' +(C2_1) e~ 4 (e 




FIG. 2: Sensitivity of the renormalised Yukawa coupling h v to 
the two-body Yukawa coupling h v fi at fixed (a) UV and 

(b) IR flow of the Yukawa coupling. We use k = fep = leV, 
abgfc = \ip,o/(&n) = 0.38 as appropriate for 6 Li and T = 0.5, 
c _1 = 1. The different curves correspond to the two-body 
value of the Yukawa coupling h% fi = 3.72-(10 3 , 10 4 , 10 5 ) (from 
top to bottom), where the last value corresponds to 6 Li while 
the first one is comparable to that of 40 K. Universality with 
respect to the value of the two-body Yukawa coupling is very 
strong even for a comparatively large (fixed) A^. 



For any possible fixed point, one has 



A,/, 



A,, 



hi = 327T =* = 1 + ^ , (66) 



4tt 



32tt 



and therefore obtains the fixed point condition for Q 

X^/hl as 



ci + (c 2 + 1)Q + 12Q 2 



0. 



In general, this quadratic equation has two distinct solu- 
tions, where the larger value of Q is infrared stable and 
corresponds to (A), while the smaller value of Q is un- 
stable and corresponds to (B). 

We finally include the effect of a nonzero negative 
chemical potential a. As soon as k 2 becomes smaller 
than a, the parameter y (|55|) rapidly goes to zero. Then, 
only the first term in the flow equations (|54[) matters. In 
this range of k, the couplings A^, h 2 ^, rh^ remain constant. 
A negative a acts as an additional infrared cutoff, such 



that the flow is effectively stopped for k 2 < —a. This 
yields a simple rough picture: the couplings \^,h 2 ,m 2 



flow according to Eq. 
(67) stops when k 2 < —a. 



for k > — g until the flow 



The fixed points are relevant not only for T = 0. We 
demonstrate the influence of the fixed point (A) on the 
flow of the Yukawa coupling for T = 0.5 and c = 1 in 
Figs. El [3] In Fig. [2] we observe that the renormalised 
Yukawa coupling h v at small k becomes almost indepen- 
dent of its initial value h 2 , if the latter is large. Figure[3] 
demonstrates the influence of the background scattering 
length cibg for a fixed value of a, again for T — 0.5 and 
c = 1. 
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FIG. 3: Sensitivity of the renormalised Yukawa coupling h v 
to the two-body background scattering length \^ for fixed 
two-body Yukawa coupling h v ,o = 3.72 • 10 5 , k — fcp = leV, 
as appropriate for 6 Li: (a) UV and (b) IR flow of the 
Yukawa coupling. We plot curves for abgfc = \^/(8n) = 
-0.38-(lCT 3 , 1CT 2 , 1CT 1 , 1) (the last value corresponds to 6 Li) 
and T = 0.5, c _1 = 1. The uppermost line in (b) has the 
largest Gtb g and starts in (a) as the lowest line. 



IX. DENSITY AND CONDENSATE FRACTION 



Let us next compute the density 



n 
fc 3 



du 



(68) 



where the a derivative will be taken at fixed p and 
k. This is possible, since for the minimum dup/da = 
du/da\ Po + (du/dip)(d(fo/da) = du/da\ Po . Before doing 
a more detailed calculation, we approximate the poten- 
tial by 



1 



u(p, a) = u(a) + rh v {a)p + -\p 



(69) 



with X v independent of a. The density then receives 
contributions from unbound atoms and molecules (hp + 
2«m) and atoms in the condensate (he) f° r /°o 7^ 0, 



hp + 2fiM 



du 

"dB 



nc 



dm 2 
da 



(70) 



For an estimate of the condensate fraction f^c = nc/n, 
we use (up to small corrections involving K- ln ) 



do 



= -2Z« 



1 + O [cqK In 



(71) 



where we have used (l4"4")l and (pjQ)) . For small K, the term 
~2Z V dominates and we conclude that the condensate 
fraction is given by the renormalised order parameter po, 



n c = 



2p 



(72) 



The above formula is one viable definition of the conden- 
sate fraction. In general, it still needs to be related to 
the definition of a condensate fraction as measured in a 
particular experiment. 

The density receives contributions from atoms and 
molecules with different momenta. In order to describe 
each momentum mode accurately, we first compute the 
change of the density with k, 



dpda 



dpo 
dk 2 ' 



d A d duo d duo 
dk 2 dk 2 da dk 2 da 

Pu • 

(73) 

and then integrate from k ln to k = 0. The r.h.s. will 
be dominated by momenta q 2 rs k 2 . The first term in 
Eq. ((721) is given by the a derivative of Eq. (J27J. The dif- 
ferent contributions have an approximate interpretation 
in terms of unbound atoms, molecules and the conden- 
sate density 



dhp 

dh^ ~ ' 

dn c 
dk 2 



d_ 

da' 
d 2 u 



dh 



dpda 



Pa 



dpo 
dk 2 



dk 2 
VvPo 



-^tCb(po), 



(74) 



For example, for po = ft one rediscovers the integral over 
the usual Fermi distribution (Q^ = 1) 



dfl F = JLr (ft) 



dk 



(75) 



d 3 q 



dk 2 
_d_ 

dk 2 J (2tt) 2 exp{(A/,<7 2 + k 2 - a)/f] + 1 ' 



In ([75]) we have traded the a derivative for a k 2 derivative, 
since 7 depends only on the combination k 2 — a. As long 



as Aj, depends only on K 



a, the k derivative 



also acts on A^, and Eq. ([75)1 can be integrated, 



up 



d 3 c 



(2tt)3 exp{{A 4 ,q 2 - a)/f} 
+h F (k in ) -£ F (p = 0, fc in ). 



(76) 



We may evaluate the initial value hp(k- ln ) in perturba- 
tion theory, finding both quantities being exponentially 
suppressed, hp(k in ) — (f(p — 0, k in ) sa 0. In our approx- 
imation, (A.jf, = 1) hp is simply the density of a free gas 
of fermionic atoms as long as po — 0. In the presence 
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of a condensate (po 7^ 0) , we typically find a flow where 
Po{k) = for k > k c , such that the flow of h-p(k) remains 
unchanged in this range. On the other hand, the contri- 
bution from modes with q 2 < k 2 will be suppressed by 
the presence of a gap A = h v ^/po in the propagator. We 
emphasise that the derivative d/dk 2 in Eq. (f75|) does not 
act on poj since the a derivative in Eq. (174j) is taken at 
fixed p. Therefore the flow for hp has to be integrated 
numerically if po 7^ 0. 

For the computation of the molecule density jim, we 
need d(s/da. Let us concentrate on p — and neglect 
the a dependence of Q v and such that 

^ - on [ d3 $ 9 
da ~ ^ J (2tt)3 da 

(exp {{A v q 2 + 2k 2 + ml)/f) - l) ~* 



dm 2 , d f d 3 q 



^ v da dk 2 J (2tt) 3 

x (exp{(A v g 2 + 2k 2 + m%)/f} - 1 

We next use Eq. ^ 

.2 d 



(77) 



dm* 
da 



-2-m v — hiZ v . 



(78) 



The second term can be approximated for small K by 







m 2 — — In Z a 



'da 



2(k 2 -a)' 



(79) 



If we neglect this term and set Q v — 1 we find the intu- 
itive formula 



n M 



(27r) 3 exp{(^<7 2 + m2)/T} -1 



(80) 



This is the expression for a free boson gas. What is cru- 
cial, however, is the appearance of the inverse propagator 
for the renormalised bosonic field ip instead of the micro- 
scopic or "bare" molecule field (p. The propagators for 
the renormalised and bare fields are related by a rela- 
tive factor Z v , such that we find for the density of bare 
molecules 



«M,ba 



= z~ 



(81) 



As a consequence, one may have a substantial molecule 
fraction 2hyi/n even if the density of bare molecules is 
tiny. These features reproduce the results from a solu- 
tion of the Schwinger-Dyson equations @ where a bare 
molecule density in accordance with the experimental re- 
sult of Partridge et at. [32J has been found. It is one of 
the important advantages of our functional renormalisa- 
tion group approach that it accounts for the distinction 
between renormalised and bare molecules in a very direct 
and straightforward manner [l4| . Within a Hamiltonian 
formulation, this issue is related to the mixing between 



open-channel and closed-channel atoms, and its correct 
treatment is crucial for a quantitatively reliable descrip- 
tion of the crossover. We emphasise that the renor- 
malised molecule field plays an important role even on 
the BCS side of the crossover. The composite bosons are 
crucial effective degrees freedom, even though the micro- 
scopic theory can be very well approximated by a point- 
like interaction between fermionic atoms without any ref- 
erence to molecules. On the BCS side, the renormalised 
field <p describes Cooper pairs. Nevertheless, one never 
needs these physical interpretations explicitly, since the 
density only involves the a dependence of the potential 
at its minimum, which is a well-defined quantity. 

We finally turn to the condensate density. If we can 
approximate the a dependence of (u—u) by a term —lap, 
as suggested by Eq. (178)) , we infer 



Ilk 2 



= 2 



dpo 
dk 2 



Zf/pPo — 2Z V 



d 



Po 



dk 2 Z u 



(82) 



(This approximation needs to hold only in the range of 
k where po differs from zero.) The dominant flow of he 
typically arises from a region where the k dependence 
of Z v is already sub- leading, resulting in he ~ 2po, a s 
found above. In practice, all these various approxima- 
tions of our analytical discussion need not to be made, 
since it is sufficient to follow the flow of h(k) numerically, 
starting from an initial value h(ki n ) « and extracting 
the physical density for k = 0. 



X. VACUUM AND TWO-ATOM SCATTERING 

It is one of the advantages of our method that it 
can access simultaneously the many-body physics of a 
gas in thermal equilibrium and the two-body physics of 
atom scattering and molecular binding. Indeed, the two- 
body physics describes excitations above the vacuum. In 
turn, the vacuum and the properties of its excitations 
obtain in our formalism simply by taking the limit of 
vanishing density and temperature. (More precisely, the 
limit should be taken such that the ratio T/h 2 ^ 3 is large 
enough that no condensate occurs.) The scattering cross 
section between two atoms can then be directly inferred 
from the Yukawa coupling and the propagator of the 
molecule field in the vacuum. 

For T = 0, the condition of zero density requires a < 
in order to ensure Uf = 0, cf. Eq. |751) . On the other 
hand, we infer from Eq. (|5D|) that m 2 ^ > is needed for 
= 0, whereas a vanishing condensate density requires 
po = 0. The vacuum is the state that is reached as density 
and temperature approach zero from above. It therefore 
corresponds to the boundary of the region where a < 
0, m 2 > 0. There are two branches of vacuum states 
(for a more detailed discussion, see [I3|)- The first has a 
negative a — a a < and m 2 = 0. In this case, the single- 
atom excitations have a gap — a a > which corresponds 
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to half the binding energy of the stable molecules, 
e M = k 2 a A /M , e M = 2<t A - 



(83) 



For the molecule phase, the effective interaction depends 
on the energy uj of the exchanged molecule 



We identify this state with the "molecule phase" where 
stable molecules exist in the vacuum. This phase is re- 
alized for B < Bq, with em or a a being a function of B 
that vanishes for B —> Bq. The other branch corresponds 
to a = 0, m 2 > 0. This "atom phase" of the vacuum is 



A^,a(w = 0,p = 0) 



Avj,r(w = -2&a,p = 0) = - 



2Z v (&a)(ta ' 
4Z ip (& A )a A ' 



(90) 



realized for B > B , with a function of B vanishing We find for the atom phage 



for B — > B . In the atom phase the "binding energy' 
£m = k 2 m 2 /{2M) is positive and the "molecules" corre- 
spond to unstable resonances. We observe a continuous 
phase transition between the molecule and atom phase 
6, 14] for m 2 = 0, a = 0, corresponding to the location 
of the Feshbach resonance at B — Bq. This fixes 8v in 
the initial value of m 2 (TT5)) by the requirement that for 
v = the mass term to 2 (cr) vanishes precisely for a = 
if T = 0. 

The vacuum state can be used to fix the parameters of 
our model by direct comparison to experimentally mea- 
sured binding energies and cross sections. Let us first 
consider the molecular binding energy cm(-B) that can 
be computed from Eq. ([50|) by requiring m 2 = 0, i.e., 



CTA(B)k 2 




(84) 



In the vicinity of the Feshbach resonance (small |cm|) the 
binding energy depends quadratically on B 



£m — 



(167r) 2 M/x 2 (B - Bf 



Using 



ao = 



8?rfc 



kh 2 

a V>,0 ~ 4M 2 ' 



(85) 



4?r p,(B -B )-e 



«o 



Eq. 



reads 



(86) 



(87) 



The scattering length for the fermionic atoms is deter- 
mined by the total cross section in the zero-momentum 
limit, a — Aira 2 , such that 

_ M- A^cff 

a — —A J, G fr — — . (SSJ 

4tt v ' 87rfc 

More details are described in the appendix. For the atom 
phase (eM > 0), one has 



m 2 (a = 0) 



A^,,o + A^,r. (89) 



ao 



Anfl(B - B Q ) 



«o 



res i 

(91) 

0. Eq. dHU) therefore 



which agrees with a (|86[) for cm 
relates the binding energy in the molecule phase to the 
scattering length in the atom phase. 

The value of h Vt Q can now be extracted from the 
resonant part, a rcs , whereas A^o follows from the B- 
independent "background scattering length" ao = ab g 



which has been measured as a 



(Li) 
bg 



-1420a B 



i (K) 



174ae (ae the Bohr radius), or, expressed in units 
h = c = k B = 1 (a B = 2.6817 • lO^eV- 1 ), agf = 
= 4.67 • 10~ 2 eV _1 . For k = leV cor- 



-0.38eV" 



i (K) 



responding to a density n — 4.4 • 10 12 cm 3 , we find the 
dimensionless expressions for 6 Li and 40 K: 



£2(Li) _ 
<p,o ~ 

A (Li) _ 



3.72- 10 5 
-9.55 , 



■2(K) 
V.O 



6.1 -10 3 , (92) 



A 



(K) _ 



1.17. 



From these values, one can compute the initial values 

fcin = 10 3 , 



A^,in) h 2 „ For 6 Li, we use the values 



aSS = - 2 - 51 • 1(r2 < 



(93) 



whereas for 40 K we take 

fc in - V300, = 1.67 • 10 5 , Xfl = 6.14. (94) 



For 40 K, we observe large values of h^ in and A^,,i n . Sim- 
ilar values for 6 Li are observed in the corresponding K 
range, e.g. h^^K 2 = 300) = 6490, xf i] (K 2 = 300) = 
2024. The broad Feshbach resonances indeed describe 
strongly coupled systems. In our approximation, the so- 
lutions of the flow equations for /i 2 , A,/, and Z v do not 
depend on m 2 provided m 2 remains positive during the 
flow. 



XI. SCATTERING LENGTH FOR MOLECULES 

So far, we have only discussed the contributions of the 
molecule fluctuations to du/da in order to extract the 
molecule density um- In this section, we discuss their 
influence on the effective potential for T = more sys- 
tematically. In particular, we will extract the scatter- 
ing length for molecule-molecule scattering in vacuum 
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from X v — d 2 u/dp 2 \po. From Eq. (|26[) . we infer for 
the derivative of the potential (primes denote derivatives 



with respect to p) for the simple case of Q v 

Qip = 1) 



1 and 



A/> — Zip 



dk 2 



dx- 



4vr 2 J y / A v x 2 +2k 2 +u , + 2pu" ^A v x 2 +2k 2 - 



pu 



' 2 



A v x 2 - 



2k 2 - 



3pu" 2 + 2p 2 u"u" 



A v x 2 - 



2k 2 



2pu" 



.(95) 



r 



This contribution vanishes for p = 0, and we find in the 
symmetric phase no influence on the running of m 2 . 
The bosonic contribution to the running of X v — u"(0), 



dk 2 V 2lT 2 ' 



{A v x 2 + 2k 2 + m 2 ) 2 (96) 



\' 2 



8^ A 3 J 2 (2k 2 + m 2 ) 1 / 2 '' 



which approaches R a in the infrared limit k — > 0, R(k — 
0) = R a . With K 2 = k 2 — a, we find the flow equation 



K" 



dR 
dK 2 



R 



R 



= R 



8K 4 X V dh% 
h% dK 2 



8K 4 dX u 



h% dK 2 



(100) 



8^ 



2V2X 2 V K 4 



3 hi 



32n K nh 2 v J 2k 2 +ml 



h 2 



, -R-3 
32irK \ 2 



k 2 -. 



k 2 + ml/2 



R z 



combines with the fermionic contribution and the anoma- 
lous dimension 



j9 h lK (;2 -x-3/2 , 9 \(F) , 9 (B) 

dk 2 v ~ 32it 1 a) + dk 2 ^ + dk 2 v ' 



(97) 



_d_ (F) _K r H * 2 

dk 2 v 8tt 2 J (x 2 + k 2 -a) 4 

o 

= _J^. {k 2 -o)^l 2 . 

Similar to the flow equation for A^, the equation for the 
gradient coefficient A v decouples from eqs. (|38I46[) in the 
symmetric phase at T — and reads 



(98) 



There are no boson contributions in this regime. The 
equation is solved by A v — 1/2 for our initial condi- 



tion A„ 



1/2 (the initial condition is actually irrele- 



vant, since A v = 1/2 is an IR attractive fixed point of 



Eq. (USD). 

It is interesting to compare the scattering length for 
molecules om = X V M / '(2tt) with the fermionic scattering 
length a (|55j) . For this purpose, we are interested in the 
flow of the ratio i? a = olm/cl — 2A y /Ai/, )e ff- For the sake 
of a simple analytic discussion, we take A v = 1/2, A^,,o = 
and consider the broad-resonance limit. Let us first 
consider 



R 



8X v {k 2 



(99) 



Inserting the fixed point behaviour h 2 = 32irK yields 



K 



dR 
dK 2 




+ m%/2 



(101) 



If we neglect the contribution from the molecule fluctu- 
ations (the last term in Eq. (|101|) ~) we find an infrared 
stable fixed point at R = 2. For this fixed point, one 
obtains for k — > 0, a — ► a a- 



X, 



4o- A 



2em 



(102) 



to be compared with A^, j0 ff = — h 2 /{2eyi), as appropri- 
ate for the molecule phase. For the atom scattering 
in the molecule phase, the propagator of the exchanged 
molecule has to be evaluated for nonzero u — —iyi/2, 
cf. Eq. (f9T)]) . We discuss this issue in the appendix. We 
conclude that this fixed point corresponds to the mean- 
field result om = 2a, R a = 2, see also [14j . 

The bosonic fluctuations will lower the infrared value 
of R. Using the fixed point for h 2 ^ and the expression 
([84]) for the binding energy in Eq. (|49[) . one finds 



-f=,k 2 



V^<W2fc 2 - e M . (103) 



At the resonance point £m = 0, the flow equation for R 
becomes 



dR 



= -3 



3 :R 



--R 2 



dk 2 2 ' V3 

The infrared-stable fixed point now occurs for 



R* 




1.325. 



(104) 



(105) 
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Precisely on the Feshbach resonance (B = Bq), the 
asymptotic behaviour obeys the scaling form 



A. 



(106) 



For B < Bq, the mass term m 2 is smaller as compared 
to the critical value for B — Bq. For a given K the 



expression y 2k 2 + to 2 , — y 2K 2 + 2a + to 2 , is therefore 

smaller and the term ~ R 2 in the flow Eq. is enhanced. 
As a consequence, \ v turns out smaller than the critical 
value ([TPS]) . 

It is instructive to study the flow in the range where 
k 2 <C — cm/2, K w yj— £m/2 + fc 2 /V~ 2e M and therefore 



fc 2 + ^ M 2fc 2_ 



The flow equation for R reads 



(107) 



d - 



3^ 



V2fc 2 - e 



2A- 



M 7? 2 



For k —> 0, the second term dominates, 



fc 2 -e M /2' 
(108) 



<9i? 



V2R 2 



dk ^F-e M /2 
and the running stops for A; — > 0, resulting in 

' V2fc tr 



R^ 1 (k = 0) w ^ _1 (fc tr ) + \/2Arsinh 



V— cm 



(109) 



(110) 



Here, fct r is a typical value for the transition from the 
scaling behaviour (|104|) for k 2 » — £m/2 to the boson 
dominated running ()109|) for k 2 <C —em/2. We may give 
a rough estimate using fctr = y— cm/2 and R(kt r ) = R*, 



R(0) 



R* 



1 + V2Arsinh(l)^ 



0.5. 



(Ill) 



The value for a numerical solution turns out somewhat 
higher and we find for k — > 



R(k^0) = ^ =0.81. 
a 



(112) 



This is mainly due to the presence of fermion fluctuations 
which enhance R. 

Following the flow of X v in a simple truncation thus 
predicts ayi/a ~ 0.81, in qualitative agreement with 
rather refined quantum mechanical computations [33j ] . 
numerical simulations [9(, and sophisticated diagram- 
matic approaches (34j |. ayi/a — 0.6. The resummation 
obtained from our flow equation is equivalent to the one 
performed in (3f|, yielding om/» — 0.78 ± 4. Functional 
RG flows in the present truncation with optimised cutoffs 
[HI, HI] lead to a M /a ~ 0.71 [3]. It should be possible 



to further improve the quantitative accuracy by extend- 
ing the truncation, for example by including an interac- 
tion of the type ip'ipcp*ip corresponding to atom-molecule 
scattering. Already at this stage, it is clear that the fi- 
nal value of i? a results from an interesting interplay be- 
tween the fermionic and bosonic fluctuations. In a model 
describing only interacting bosonic particles (atoms or 



molecules) , the flow of A„ can be taken from Eq. 



irr 



0,A U 



~~dk~ 



— \ 2 

2vA 



for 



(113) 



(For an application to bosonic atoms, one should recall 
that our normalisation of the particles corresponds to 
atom number 2 and mass 2M.) The bosonic fluctuation 
effects have the tendency to drive \ v towards zero, ac- 
cording to the solution of Eq. (|113[) for k — ► 0, 



A 



1 + <\/?,in&in/27T 



(114) 



This solution reflects the infrared-stable fixed point for 



A^ — \^pk 7 



resulting from 



(A v )* — 0, 



k d \ — A + A ' 
k dk K ~ K+ 2^- 



(115) 



(116) 



On the other hand, the fermion fluctuations have the 
opposite tendency: they generate a non-vanishing X v 
even for an initial value = and slow down the 



decrease of X v as compared to Eq. (|113p in the range of 
small k. The resulting R a results from a type of balance 
between the opposite driving forces. If the scaling be- 
haviour (|104p is approximately reached during the flow, 
the final value of i? a is a universal ratio which does not 
depend on the microscopic details. This is a direct conse- 
quence of the "loss of microscopic memory" for the fixed 
point (fTUS) . 



XII. UNIVERSALITY 

The essential ingredient for the universality of the 
BEC-BCS crossover for a broad Feshbach resonance is 
the fixed point in the renormalisation flow for large /i 2 
(fixed point (A) in sect. IIX[) . This fixed point is infrared 
stable for all couplings except one relevant parameter cor- 
responding to the detuning B — Bq. For T = 0, this fixed 
point always dominates the flow at B = Bq as long as ef- 
fects from a nonzero density can be neglected. However, 
for any the flow will finally deviate from the scal- 

ing form due to the occurrence of a condensate po ^ 0. 
Typically, this happens once k reaches the gap for single 
fermionic atoms, k « A = h v ^/po~. (If we use k = k-p the 
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relevant scale is k = 16f2 c /(37r).) For k <C A, the contri- 
bution from fermionic fluctuations to the flow gets sup- 
pressed. Also the contribution from the "radial mode" 
of the bosonic fluctuations will be suppressed due to a 
mass-like term m^, ss 2\ v p$. Only the Goldstone mode 
corresponding to the phase of ip will have un-suppressed 
fluctuations. 

For nonzero density it is convenient to express all quan- 
tities in units of the Fermi momentum k-p and the Fermi 
energy ep, i.e., to set k = k F . As compared to the scaling 
form at zero density the flow equations are now modi- 
fied by two effects. They concern the influence of po ^ 
and the nonzero value for a — a* which corresponds to 
B = Bq. (Typically, this value <r* is positive, as op- 
posed to <t* = for the zero density case.) Nevertheless, 
all these modifications only concern the flow for small 
k. If the flow for large k has already approached the 
fixed point close enough no memory is left from the mi- 
croscopic details. We can immediately conclude that all 
physical quantities are universal for B = Bq. 

The same type of arguments applies for nonzero tem- 
perature. For T 7^ 0, a new effective infrared cutoff 
is introduced for the fermionic fluctuations. Also the 
contributions from the bosonic fluctuations get modified. 
Again, this only concerns the flow for k <C ttT. The loss 
of microscopic memory due to the attraction of the flow 
for large k towards the fixed point remains effective. 

In a similar way, we can consider the flow for B ^ Bq. 
Away from the exact location of the resonance, the ob- 
servable quantities will now depend on the relevant pa- 
rameter. The latter can be identified with the concen- 
tration c = ak F @ ■ Still, the deviation from the scaling 
flow only concerns small k, namely the range when the 
first term in Eq. (|49p becomes important. This happens 
for 



k 



16irMp,\B - B Q 



1 



h 2 k 2 C 



(117) 



If the flow for large k has been attracted close enough to 
the fixed point all physical quantities at nonzero density 
are universal functions of two parameters, namely c and 
T = T/ep. These functions describe all broad resonances. 



XIII. DEVIATIONS FROM UNIVERSALITY 

The functional flow equations also allow for simple 
quantitative estimates for the deviations from universal- 
ity for a given atom gas. Typically, these corrections 
are power suppressed ~ (Air,/Atjv) p - The microscopic 
scale Atjv roughly denotes the inverse of the character- 
istic range of van der Waals forces and we may associate 
Atjv — 1/(100<2a) with a a the "radius" of the atoms. 
More precisely, Auv corresponds to the scale at which the 
attraction to the fixed point sets in. In some cases, this 
may be substantially below the van der Waals scale, as for 
the case of 40 K, as we will argue below. The relevant in- 
frared scale Air depends on the density, Air = L(c, T)fcp- 



In view of the discussion in the preceding section, we ap- 
proximate L by the highest value of k where the deviation 
from the scaling form of the flow equations sets in, i.e., 



L = max{ | c 



(118) 



We note that a and A = h v ^/po depend on c and T. 

The power p of the suppression factor reflects the 
"strength of attraction" of the fixed point. More pre- 
cisely, we may consider the flow of various dimensionless 
couplings on . If we denote their values at the fixed point 
by aj*, the flow of small deviations Sex, — on — is 
governed by the "stability matrix" Aij 



(119) 



The relevant parameter corresponds to the negative 
eigenvalue of A. All other eigenvalues of A are positive, 
and p is given by the smallest positive eigenvalue of A. 

Typical couplings on are (A,/,, fit, mL, R). For these 
couplings, the stability matrices reads for fixed point (A): 



.4 




(120) 



with ei = (R\ - 6 + 2i?2)/(327r), e 2 = -Rl/(6y/$), e 3 = 
3 + 4R*/V3- The eigenvalues are (1,1 — l,e 3 ) and we 
conclude p = 1 . This situation is not expected to change 
if more couplings like A v , A^ or the p 3 term in the po- 
tential are included. 

We may use our findings for a rough estimate for the 
range in c" 1 and T for which deviations from universality 
are smaller than 1%. From our previous considerations 
this holds for 



< 



Aw 
lOOfcp 



A uv 
300fc F 



(121) 



As a condition for Auv, we require that all couplings 
are in a wider sense in the "vicinity" of the fixed point 
at the scale fcuv = Auv/^f- Of course, one necessary 
condition is Ayy < A. Also, the flow equations must be 
meaningful for k < k\jy- While this poses no restriction 
on Li where A^o < 0, the case of K with positive A^o 
requires fcyvco < 1 or fcjjv 20. If we would try to 
extrapolate the flow to even higher k the value of h^k) 
would diverge as can be seen from Eq. (|4Tj) . As a further 
condition for the flow to be governed by the universal 
behaviour of the fixed point we require Z v — 1 > 1. This 
tells us that the notion of a di-atom state <p starts to be 
independent of the detailed microscopic properties. For 
Li and K, this happens for scales only somewhat below 
k\ n . 

A reasonable estimate at this stage may be kyy = 

800, fcyy = 15. One percent agreement with the univer- 
sal behaviour would then hold for 



Li : | c— 1 1 < 8 , VT < 3, 
K: Ic- 1 ! < 0.15 , ^^0.05, 



(122) 
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which corresponds to 



Li 

K : 



|afc F | > 0.13 
\ak F \ > 6.7 , 



T/T F < 9, 
T/T F < 0.0025. 



(123) 



Already at this stage, we conclude that the universal be- 
haviour for K covers only a much smaller range in c and 
T as compared to Li. In particular, experiments with Li 
are indeed performed within the universal regime which 
extends far off resonance and to temperatures well above 
the quantum degenerate regime T/T F rs 1, while for 
K it will be hard to reach the truly universal domain, 
since the lowest available temperatures range down to 
T/T F ~ 0.04, and the magnetic field tuning is too low to 
resolve a regime where \ak F \ > 6.7. 

At first sight, the limitation to a rather small value 
A'liv = 15 for 40 K may seem to be of a technical nature, 

enforced by the breakdown of the flow for Aw, for k > fc^y . 
However, this technical shortcoming most likely reveals 
the existence of some additional physical scale close to 



„(K) 

V UV ' 



as, for example, associated to a further nearby res- 
onance state not accounted for in our simple model. In- 
deed, if our model (without modifications and additional 
effective degrees of freedom) would be valid for k 3> k^J , 
one could infer an effective upper bound for A,/, from its 
flow. Starting at some fc; n with an arbitrarily large pos- 
itive A^in, the value of \^(k) would be renormalised to 
a finite value, bounded by Ay (A = 0) < 87r/A"i n (cf. 
Eq. (|39|) ). This would lead to a contradiction with the 
observed value unless K m is sufficiently low. The ob- 
served value of A,/,^ therefore implies the existence of a 
scale near k^J where our simple description needs to be 
extended. This issue is very similar to the "triviality 
bound" in the standard model for the electroweak inter- 
actions in particle physics. 

As further concern, one may ask if ft 2 (fcuv) is al- 
ready close enough to the fixed point value 327r and if 
A^(fcuv) is close enough to zero. This is an issue, since 
we know of the existence of a different fixed point for 



narrow Feshbach resonances at h~„ 



A,; 



0. At the 



"f. ~ * . 

scale k\jy , the flow should not be in the vicinity of this 
"narrow-resonance fixed point" anymore, but at least 
in the crossover region towards the "broad-resonance 
fixed point" . If we evaluate the couplings at the scales 
fc£$ = 800 , 4v = 15 we obtain for cT 1 = 



Li : ftj(fcuv) - 5.00 ■ 10- 3 , A^(fcuv) 

m£(fcuv)=6.07-l(n a ; 
K : h 2 v {kuv) = 4.53- 10 3 , Ay(fc uv ) = 58.8, 

m 2 (fe uv ) = 54.0. 



-25.0, 



(124) 



For K, the relevance of the broad-resonance fixed point 
seems plausible, but the small value of /i 2 (fcuv) f° r Li 
may shed doubts. However, one should keep in mind 
that for a strong enough Ay, jC ff the distribution between 
Ay, and — /i 2 /m 2 concerns mainly the dependence of the 



scattering length on the magnetic field B and not so much 
the physics for a given a or c. Indeed, we could absorb 
Ay by partial bosonisation (Hubbard-Stratonovich trans- 
formation) in favour of a change of /i 2 and m 2 . Up to 
terms ~ g 4 , which are sub- leading for the low- momentum 
physics, a model with given Ay, ft, 2 , m 2 is equivalent to a 

model with Ay, = 0, but new values h£ and m^ 2 related 
to the original parameters by 



hi = hi 



.'2 



2Aym 2 

Aym 4 
h% " 



V'V 

K 



(125) 



The value of the new Yukawa coupling for Li, ft^fcuv) = 
4.65 -10 2 , is much larger and seems acceptably close to the 
broad resonance fixed point. Using partial bosonisation 
and the technique of rebosonisation during the flow [3lj ] , 
it may actually be possible to treat Ay, as a redundant 
parameter, thus enlarging the "space of attraction" of the 
fixed point (A). In this context, we note that we should 
include the contribution of the bosonic fluctuations to 
the running of Ay for our truncation. This will shift the 
precise location of the fixed points. In a language with 
rebosonisation where Ay, remains zero, these additional 
contributions will be shifted into the flow of m 2 and ft. 2 . 

Our estimate for the deviations from universality con- 
cerns only the overall suppression factor. Detailed pro- 
portionality coefficients co for the deviation of a given 
observable O from the universal strong-resonance limit, 
SO = cqAir/Auv, depend on the particular observ- 
able. The deviations from universality should therefore 
be checked by explicit solutions of the flow equation for 
Li and K. 



XIV. CONCLUSIONS 

In this paper, we have performed a functional renor- 
malisation group study for ultracold gases of fermionic 
atoms. We have concentrated on four parameters: the 
Yukawa coupling of the molecules to atoms h v , the off- 
set between molecular and open channel energy levels 
m 2 which is related to the detuning, the background 
atom-atom-interaction Ay, and the molecule-molecule in- 
teraction A^,. This system exhibits a fixed point for the 
rescaled couplings which is infrared stable except for one 
relevant parameter. This parameter can be associated 
with the inverse concentration c _1 = (a/cp) -1 . 

Whenever for a given system the trajectories of the 
flow in coupling-constant space approach this fixed point 
close enough, the macroscopic physics looses the mem- 
ory of the microphysical details except for one parame- 
ter, namely c _1 . In consequence, for a certain range in 
c _1 around and for a certain range in temperature T 
all dimensionless macroscopic quantities can only depend 
on c and T. Here dimensionless quantities are typically 
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obtained by multiplication with appropriate powers of k F 
or ep. 

The macroscopic quantities include all thermodynamic 
variables and, more generally, all quantities that can be 
expressed in terms of n-point functions for renormalised 
fields at low momenta. In particular, the correlation 
functions for atoms and "molecules" as well as their in- 
teractions depend only on c and T, independently of the 
concrete microscopic realization of a broad Feshbach res- 
onance. Here "molecules" refers to bosonic quasi parti- 
cles as collective di-atom states which may be quite differ- 
ent from the notion of microscopic or "bare" molecules. 
In a certain sense this situation has an analogue in the 
universal critical behaviour near a second order phase 
transition. In contrast to this, however, the universal 
description includes now a temperature range between 
T = and substantially above the critical temperature, 
and the universal quantities depend on an additional rele- 
vant parameter c _1 . We may also compare to a quantum 
phase transition at T = which is realized in function of 
a parameter analogous to c _1 . In our case, universality 
is extended to nonzero temperature as well. 

The "broad resonance fixed point" is not the only fixed 
point of the system. Another "narrow resonance fixed 
point" at /i 2 , = 0, \ v = allows for an exact solution 
of the crossover problem for narrow Feshbach resonances 
0, [36[ and a perturbative expansion around this solution. 
The narrow-resonance fixed point has two relevant pa- 
rameters c _1 and /i 2 . A typical flow away from this fixed 
point at small enough c _1 and A,/, describes a crossover 
towards the broad resonance fixed point. 

For a given physical system characterised by some 
microscopic Hamiltonian, it is important to determine 
how close it is to the universal behaviour of the broad- 
resonance limit. We have performed here a first estimate 
for the range in c _1 and T for which universality holds 
within one percent accuracy for the experimentally stud- 
ied Feshbach resonances in 6 Li and 40 K. It will be an 
experimental challenge to verify of falsify these predic- 
tions of universality. 



APPENDIX: SCATTERING LENGTH AND 
CROSS SECTION 

In this appendix, we collect useful formulae for the 
scattering of molecules in vacuum. In quantum field the- 
ory, the scattering cross section for identical nonrclativis- 
tic bosons is given by 



(A.l) 



The factor 1/2 is a convention, motivated by integra- 
tion over half the spatial angle in order to avoid double 
counting of the identical particles. The differential cross 
section is related to the scattering amplitude by 



da 



16tt 2 



(A.2) 



such that, in case of isotropic (e.g. s-wave) scattering 



Stt 



(A.3) 



Next, we relate the above results to nonrelativistic quan- 
tum mechanics: The scattering length a is a quantity 
which is directly meaningful in nonrelativistic quantum 
mechanics only. It measures the strength of the 1 jr de- 
cay of the scattered wave function at low energies as a/r. 
Equivalently, a is defined by the I = phase shift for the 
scattered wave function. This definition leads to the re- 
lation between scattering length and cross section (iden- 
tical bosons) 



o"b = 87raj 



(A.4) 



We can use Eqs. (|A.3p . (|A.4p in order to relate the scat- 
tering amplitude to the scattering length or rather to 
define this relation, 



87ra B 2A V 



(A.5) 



The scattering amplitude is given by the effective four- 
boson vertex, as obtained by functional derivatives of the 
effective action, Z^Ab = — 2\ v . In the limit of 
a point-like interaction, A^, is a constant. (Note that 
the notion of scattering amplitude is used to describe 
different quantities in quantum mechanics and QFT.) 
For fermions, the definitions are similar, 



Of 



dn- 



da 
> 

' 'dn' 



(A.6) 



Now, the integration covers the full space angle, since the 
fermions can be distinguished by their spin. The differ- 
ential cross section is related to the scattering amplitude 
by 



da _ M<*\A F \ 2 
dn ~ 16tt 2 ' 

such that for isotropic scattering 

M§\A F \ 2 



a F 



4tt 



(A.7) 



(A. 



For distinguishable fermions, one has in quantum me- 
chanics 



a = Aira 



F: 



(A.9) 



such that the scattering amplitude and scattering length 
are related by 



4ira F - . . 



(A.10) 



The scattering amplitude A F is given by the tree level 
graph for the molecule exchange fermions plus a contri- 
bution from the fermionic 1PI vertex A^. The energy- 
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and momentum-dependent resonant four-fermion vertex 
generated by the molecule exchange reads 



(A.ll) 



-us - 



1" 

4M 



(AP v (w,£)-AP v (0,0)) 



with 



4~ 



-ui — 2a a 



4M 



(A.12) 



With q\,q2 the momenta of the scattered atoms, one 
finds for the momentum and energy of the exchanged 
molecule q = q\ + q 2 , uj = g?/(2M) + <z|/(2M) - 2a A - 
Here, we take into account the binding energy - the en- 
ergy of an incoming atom is — oa+^I /(2M). We now con- 
sider the limit qi — > and work in the broad-resonance 



limit hz, — > oo. The vacuum in the molecule phase is 



characterised by m 2 — 0, and we end up with 



ty.eff - 



A^,o — 



47T 



(A.13) 



One infers for the scattering length of atoms in the 
molecule phase 



MX^ eS (-2a A ,0) 



1 



+ a . (A.14) 



4tt V-2A/cta 
The relation between our rescaled quantities and the 



scattering amplitudes is given by 



\, = 2Mk F X v , X v = — , A^, i0ff = 2Mk F X^ te s, 



Vo = aA - ~> m <p = , ■ (A.15) 



With Mb = 2M, the ratio between molecular and atomic 
scattering length (om = ob, «m^f = A ¥ ,/(47r)) therefore 
reads 



2A V /Z 2 



2K 



a X^,, cS {uj,q) A^ !Cff (-2CTA, q) 



(A.16) 



Using the definition of the renormaliscd Yukawa coupling 



h 2 = h?/Z v , we may express 



A^, e ff(— 2cta, 0) — A^,,o 



8tt 



y/— &A 4<TA 2e M 

(A.17) 

In the last two expressions, we have used the leading 
terms in Eqs. gj), (05) for = V^oa 0. 
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